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Abstract. We study experimentally systems of orthogonal polynomials with 
respect to self-similar measures. When the support of the measure is a Can- 
tor set, we observe some interesting properties of the polynomials, both on 
the Cantor set and in the gaps of the Cantor set. We introduce an effective 
method to visualize the graph of a function on a Cantor set. We suggest a new 
perspective, based on the theory of dynamical systems, for studying families 
Pn{x) of orthogonal functions as functions of n for fixed values of x. 



1. Introduction 



The classical theory of orthogonal polynomials [Gautschi 2QQ4j [Szego 1975] al- 



lows you to start the Gram- Schmidt process with virtually any measure. In the fam- 
ilies of polynomials that are usually studied, the measure is either absolutely contin- 
uous or discrete, but the general theory allows one to use a singular continuous mea- 
sure. In recent years there has been an interest in the case of fractal measures. See 
for example [Barnsley et al. 1983^ - [Barnsley et al. 1985 ,[Bessis and Moussa 1983 , 
jMantica 1996j - [Mantica 2QQQj . In particular, G. Mantica has developed algorithms 
to efficiently compute the coefficients of the 3-term recursion relation, and hence 
the polynomials, in the case of a self-similar measure [Mantica 1996 , M antica^^QQO] . 
In this work we use these algorithms as a tool to study the polynomials experimen- 
tally, looking for interesting patterns and conjectures. We view this as part of a 
general program to explore portions of classical analysis related to fractal measures. 

Let /i be a measure on the line, and for simplicity assume that /i is a probability 
measure supported on the unit interval. The related system {Pn{x)} of orthogonal 
polynomials is characterized as follows: 

Pn{x) is a polynomial of degree n with a non- negative coefficient of x"^, (1) 
P„(.)P„(.)<iM«) = ,Ortho„o™*y. (2) 
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The general theory imphes that there is a 3-term recursion relation 

{x) + (3) 

where r^, An > 0, n G Z+ are determined by the measure /i and P_i =0, Pq = 1- 
The coefficients r^, An are also called the entries of the Jacobi matrix , an infinite 
symmetric tri-diagonal matrix J such that Ji^i = Ai and Ji,i+i = Ji+i,i = r^+i for 
2 G Z^o. Note that Eq. ^ allows us to find the polynomials recursively: 

Pn+l{x) = -^{{X - An)Pn{x)) - ^Pn-l{x) (4) 

We should point out that while in principle Eq. Q allows computation of the 
coefficients of P^, this computation could be unstable. In the case that An is 
constant (as we see below when our measure is supported on a Cantor set), it is 
often easier to compute the coefficients of P^ as a polynomials of {x — An)- When 
the An vary, we alternatively use Eq. Q to compute the values Pn{x) for specific 
X- values. 

A measure fi on the line is said to be self-similar if there exists an iterated 
function system (IFS) of contractive similarities {Pi}^i and a set of probability 
weights {pi}f^i such that 

N 

any measurable set A, or equivalently, (5) 

i=l 

j f dji = J f ^ Fi djji^ for any continuous function /. (6) 

In this paper we restrict our attention to the family of IFS's with N = 2 and such 
that 

F,{x) = y, F2(x) = + l (7) 

where > 2 is a parameter. When R = 2 and = P2 = ^ we obtain Lebesgue 
measure on [0,1], and the corresponding polynomials are essentially the classical 
Legendre polynomials. (Actually, the classical Legendre polynomials are orthogonal 
with respect to Lebesgue measure on [—1, 1] and are normalized differently, but the 
differences just involve rescaling the axes.) When R = 2 and pi 7^ P2, we refer 
to the corresponding polynomials as Weighted Legendre Polynomials (WLP). The 
measure /i is singular but not supported on any proper closed subset of [0, 1]. When 
R> 2 we will always take pi = P2 = | , and we call the corresponding polynomials 
Cantor Legendre Polynomials (CLP). The measure ja is then supported on a Cantor 
set Cr characterized by 

Cr = FiCr[JF2Cr (8) 
The standard Cantor set and Cantor measure correspond to R = 3. We refer to 
the intervals in [0, 1] \ Cr as gaps. The largest gap is the interval (;^, 1 — ;^), and 
there are 2^ gaps of length (1 — ;|). 

The behavior of Pn{x) in the CLP case is quite different on the gaps and on 
the Cantor set Cr. In order to visualize the graphs of Pn{x) on Cr we introduce 
the distorted Cantor set Cr^c (also known as the Smith- Volterra-Cantor set or fat 
Cantor set), obtained by reducing the size of the gaps by a factor of e (a parameter 
that we choose). Note that CR^e is still a topological Cantor set, but it has positive 
Lebesgue measure. There is an obvious one-to-one correspondence between Cr and 
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Cn^e that identifies regions between corresponding gaps. We use this identification 
to graph functions defined on Cr against Cr^c- 

In section |2] we present data for the entries of the Jacobi matrix. In the CLP 
case we note the different behaviors of for even and odd n, and the occurrence 
of smah values. In section |3] we display graphs of the polynomials. In the CLP case 
we show graphs of the restrictions to the Cantor set and to the gaps. 

We then discuss various features of the data. In section [4] we discuss the as- 
sociated Dirichlet kernels. By using the Christoffel-Darboux formula we are able 
to relate approximate identity behavior with small values of r^. In section [5] we 
discuss some approximate equalities relating CLP restricted to the Cantor set. In 
particular, P2n(^) is approximately equal to P2n+i(^) on the right half of Cr^ for 
large n. We also define a "shuffle" map that approximately preserves P2n+i when 
n is a power of 2. In section [g] we discuss the behavior of Pn{x) on the gaps in 
the CLP case. On the central gap, for high n, Pn{x) vaguely approximates either a 
Gaussian or the derivative of a Gaussian, depending on n mod 2. More precisely, 
we find that P2n{x) = C2ne~^^'^^"^^^ where a{x) = a2n{x) = 2 at x = 1/2, and 
a{x) ^ 2 for X around 1/2. On the other gaps, for large enough n, the behavior 
is roughly the same. In section [t] we discuss the behavior of Pn{x) at the points 
X = 0, 1/4, 1/2, 3/4, 1 in the WLP case, and we contrast the results with the known 
behavior for Legendre polynomials. This study leads to the dynamical systems 
perspective. 

Instead of thinking of Pn{x) as a function of x for fixed n, we look at Pn{x) as 
a function of n, for a fixed x. Because of the 3-term recursion relation, it is more 
natural to look at vectors (Pn(^), Pn+i(^)) in the plane. Then there exist 2x2 
matrices Mn{x) such that 

( "i^ ) = ( ) ■ (») 

In fact. 



So, we are really looking at orbits of a time-dependent linear dynamical system. 
We can then define the analogy of the Mandelbrot set as all x for which the orbit 
is bounded. It appears that this coincides, more or less, with the support of the 
measure. Then for each x in this Mandelbrot set, we can define a Julia set J{x) as 
the limit set of the orbit. We display some examples of these Julia sets. We then 
extend the investigation of section [7| to generic x for WLP and CLP in section [S] 
We end with a short concluding discussion in Section |9] 

This paper should be viewed in the context of a long term effort to under- 
stand topics in classical analysis extended to fractal measures. The following ref- 
erences are just a sampling of this work: [Bird et al. 2006) IColetta et al. 2004} 
Dutkay and Jorgensen 2006, Huang and Strichartz 2001, Jorgensen and Pedersen 2000| 



Kigami 2001, Laba and Wang 2002, .Lau and Wang 1993, Lund et al. 1998j 
[Strichartz 1990 - [Strichartz 2006 . More data may be found at 
[www . math . Cornell . edu/~orthopoly , 

Acknowledgement: We are grateful to Giorgio Mantica for allowing us to use 



his codes. 
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2. Entries of the Jacobi Matrix 



The coefficients and in the 3-term recursion relation (Eq. ([3|) determine 
the polynomials in a rather subtle way. More work is needed to clarify this relation- 
ship. In this section we report data for our two classes of examples. In Figs. [l][4]we 
graph An and versus n, for several choices of pi in the WLP case. For classical 
Legendre polynomials we have An — 1/2 for all n by symmetry about x = 1/2, and 
lim^^oo = 1/4. The WLP case shows a small but significant difference from this 
model case. 



A vs. index for WLP P =.6 




1000 2000 30OO 4000 5000 SOOO 7000 SOOO 9000 10000 

Index n 



Figure 1. Plot of An vs. index n, 1 < n < 10000 for pi = .6 



A vs. index for WLP P =.7 



1000 2000 3000 4000 5000 BOOO 7000 SOOO 9000 10000 

Index n 



Figure 2. Plot of An vs. index n, 1 < n < 10000 for pi = .7 

In the CLP case, all = 1/2 by symmetry about x = 1/2. In Fig. [sjwe graph 
Tn versus n for R = S. 

Problem 2.1: What is the nature of the limit set 



Do the measures 



n?^=i d(U~^{r„})? 

^ m 
— h ^ — ^ 



m — k 



(11) 

(12) 



n=/c+l 
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r vs. index for WLP P =.6 




Figure 3. Plot of rn vs. index n, 1 < n < 10000 for pi = .6 

r vs. index for WLP P =.7 
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Index n 



Figure 4. Plot of vs. index n, 1 < n < 10000 for pi = .7 



converge weakly to some fractal measure as m^k tend to infinity in some specific 
manner? 

A striking feature of the data is the different behavior of for n even and n 
odd. In Figs. |6]and[7|we show the same data as Fig. |5j separating the even and 
odd values of n. Another striking observation is that some values of for n even 
are close to zero. 

Problem 2.2 : What is 

liminfrn? (13) 

In particular, is the liminf zero? What is the sequence of n's along which the liminf 
is attained? 

As we will see later, having values of close to zero has interesting implications. 
We could also ask for the limsup, but it is not clear what significance this has. 
Conjecture 2.3: We always have 



1 

n _ 2 



(14) 
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r vs. index for CLP R=8 





1000 2000 30OO 



7000 3OO0 90OO 10000 



Figure 5. Plot of vs. index n, 1 < n < 10000 for CLP, R 

r (n odd) vs. index for CLP R=8 
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Index n (odd) 



Figure 6. Plot of Vn vs. index n, 1 < n < 10000 for n odd, CLP, 

r (n even) vs. index for CLP R=8 
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Index n (even) 



Figure 7. Plot of Vn vs. index n, 1 < n < 1000 for n even, CLP, 
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3. Graphs of Polynomials 

In Fig. [sjwe show the graphs of Pn{x) for 1 < n < 5 for the WLP with pi = 1/2, 
so these are just (rescaled versions of) the classical Legendre Polynomials, computed 
using G. Mantica's algorithm. In Figs. [9] and 10 we show the same functions for 
pi = .6 and pi = .7. Already we observe that the symmetry is broken in a decisive 
fashion, as these functions are much larger near x = 1 (where the measure "has less 
weight" ) than near x = (where the measure "has more weight" ) . This observation 
is expected since the Gram-Schmidt process with respect to /i almost immediately 
implies the following property: Pn{x) is the unique n^^ degree polynomial, with 
highest degree coefficient d^, which minimizes the i^^([0, l],/i) norm. (In the case 
Pi = 1/2, we have dn = YYi^i Y -'-Z^* ~ ^^16 — 4/z^.) Therefore, Pn is expected 
to be the smallest where /i "has more support." In Fig[TT]we show the graphs of 
Pn{x) for 49 < n < 52 for pi = .7. The same graphs are shown in Fig. 12 with the 
y-ax.is truncated to more clearly display the structure of the functions. 

Conjecture 3.1: For fixed pi > 1/2 and any given e > the WLP are uniformly 
bounded on the interval [0, 1 — e]. 



(Traditional) Legendre Polynomials R=2 and 50-50 Weighting 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 8. Plot of Pn{x) on the unit interval with pi = P2 = .5. 
This gives the classical Legendre Polynomials (up to renormaliza- 
tion). 



Next we look at the CLP case, where certain features of the Pn{x) increase in 

we show the graphs of Pn{x) for 1 



complexity. In Figs. 13 



and 



14 



< n < 5 on 

the whole interval for two different choices of R. In Fig. [15] we show the graph of 
^52 (^) on the whole interval for R = A. Note that the values on the central gap 
are so large that no information about the graph on the complement of the central 
gap is discernable using linear scaling in the 7/-axis. This is typical of Pn{x) for 
large even values of n = 2k and any R > 2 (though the sign of Pn{x) changes sign 
between n = 4j + 2 and n = Aj). To view these polynomials more effectively, we 
scale the y-axis logarithmically, but also multiply by the sign of Pn{x). That is, 
we plot sign(P^(x)) log(|Pn(^)|) vs x. In Fig. 16 we show the graph of P5i(x) for 



= 4, which is typical of Pn{x) for n odd and large. With the same logarithmic 
scaling. Fig. 17 shows P52(^) for = 4 restricted to the gap [1/i?^,] = [1/16, 3/16]. 



This behavior is typical of Pn{x) for large n. Note that the maximum value is quite 



STEVEN M. HEILMAN\ PHILIP OWRUTSKY^, AND ROBERT S. STRICHARTZ^ 



Weighted Legendre Polynomials with R=2 60-40 weighting 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 



Figure 9. Plot of Pn{x) on the unit interval withpi = .6,p2 = -4. 
Recall that the 60 percent weighting is given to the map which 
contracts towards zero. 



Weighted Legendre Polynomials with R=2 70-30 weighting 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 10. Plot of Pn{x) on the unit interval with pi = .7,^2 = -3. 



P^(x) forWLP p=.7 



100- 
CL* 50 



100- 
Cl" 50 



0.2 0.4 



100 
Q.^ 50 



150 ■ 
g 100 
Q.^ 50 

— - 
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-50 L 



Figure 11. Plot of Pn{x) on the unit interval with pi = .7, j92 = -3 
for 49 < n < 52. 
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P^(x) for WLP p^=,7 ^rth y-axis restricted 




0,2 0.4 0.6 0.3 1 0.2 0.4 0.6 O.S 



Figure 12. Plot of Pn{x) on the unit interval withpi = .7,^2 = -3 
for 49 < n < 52 with the y-ax.is restricted to the range —10<y< 
10. 



large, but still it is small relative to the maximum value on the central gap. We 
will discuss this behavior more in Section |6l 



P (X) CLP R=3 




0.1 0.2 0.3 



Figure 13. Plot of first five Cantor Legendre Polynomials, R = 3. 



Next, we graph CLP on the Cantor set Cr using a distorted Cantor set CR^e 



for the X-axis. In Figs. 18 and 19 we display the same functions as in Figs. 13 and 
14 Already we see that the values of Pn{x) are considerably smaller on the Cantor 
set. In Fig. |20]we show Pn{x) restricted to the Cantor set for 97 < n < 100 and 
R = 8. We will comment in detail about some of the structure of those graphs in 
Section [5] Many more examples may be viewed at [Owrutsky 2005| . 
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P (X) CLP R=3 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



X 

Figure 14. Plot of first five Cantor Legendre Polynomials, R = 
4. As in the Cantor Legendre Case, the polynomials have lowest 
absolute value where the measure has the most weight (as ex- 
pected). 



P (x) CLP R=4, logarithmically scaled 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 



Figure 15. Plot of P52(^) for CLP R = 4, logarithmically scaled 
on the entire interval. Notice the large values obtained by the 
polynomials off the Cantor set (10^^) compared to the behavior on 
the Cantor set (for I < j < n large, \Pn\ < 20). The overall shape 
in the gap is Gaussian, and this occurs for all n = Ak^ k G Z+. For 
n = 4/c + 2, the behavior in the gap is roughly a negative Gaussian. 
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Pg^(x) CLP R=4 

40 I 1 , , , 1 1 1 , , 1 




-40 I 1 ' ' ' 1 1 1 ' ' 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 

Figure 16. Plot of P^iix) for CLP R = 4, logarithmically scaled 
on the entire interval. This is roughly the derivative of a Gaussian 
on the center gap. For all n = 4A: + 3 we observe this shape, and 
for n = 4/c + 1, the negative of this shape. 



P (X) CLP R=4 on Secondary Gap 




0.06 0.08 0.1 0.12 0.14 0.16 0.18 

X 



Figure 17. Plot of P52(^) for CLP R = 4, logarithmically scaled 
on the secondary gap [1/R^{R— l)/R] = [1/16,3/16]. Again, we 
observe a Gaussian shape. 
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Cantor Legendre Polynomials on traditional Cantor Set (R=3) 



-1.5^ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 8 0.9 1 
X (Distorted Cantor) 



Figure 18. Plot of first five Cantor Legendre Polynomials, R = 3. 



Cantor Legendre R=^ 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
X (Distorted Cantor R=4) 



Figure 19. Plot of first five Cantor Legendre Polynomials, R = A 

P^(x) for CLP R=8 ^-^^^ y.^^jg restricted 




0.2 0.4 0.6 0.8 1 

X, Distorted Cantor set 



0.2 0,4 0.6 0.6 

X, Distorted Cantor set 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.6 

X, Distorted Cantor set x, Distorted Cantor set 



Figure 20. Plot of some larger n Cantor Legendre Polynomials, 
= 8, 97 < n < 100. 
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4. DiRiCHLET Kernels 



For a general function / in L^((i/i), we can expand it as a series 

CO 

n=0 

where the coefficients are given by 

= J f{y)Pn{y) dfi{y). 

The partial sums 

AT 
n=0 

may be represented as an integral 

J DN{x,y)f{y) d^{y) 
where the Dirichlet kernel is given by 

N 

Dn{x,y) = Y.Pn{x)Pn{y) 



(15) 

(16) 
(17) 
(18) 

(19) 



n=0 



The partial sums in Eq. (17) converge to / in L norm, but to get better convergence 



we need to know more about the Dirichlet kernel. In view of related results in 
[Strichartz 1993-21 IStrichartz 1994] . one might hope that there exists a sequence 
{Nk} along which the partial sums converge uniformly if / is continuous. This 
would follow by standard approximate identity arguments if we could show that 
the quantity 



\DNk{x,y)\ dfi{y) 



is uniformly bounded, and that 



lim 



/ 

J\y-x\>€ 



\DNk{x,y)\ dij{y) = 



(20) 



(21) 



for all e > 0. While we have no insight on how to establish Eq. (20), we can say 



something about Eq. (21), thanks to the Christoffel-Darboux formula 

DN{x,y) = ^-^^{PN+i{x)PN{y) - PN{x)PN+i{y)) , for X 7^ 



(22) 



y 

(It is easy to derive Eq. J22| by multiplying Eq. ([l9| by {x — y) and using the 3-term 
recursion relation Eq. (|3|.) Assuming that the polynomials Pn{x) are uniformly 
bounded on the support of /i, if we could find a sequence Nj^ such that vn^+i ~^ 0^ 
then Eq. (I^Tl) follows from Eq. (I22|). 



It appears from our data that for CLP polynomials there exist indices n such 
that Tn is close to zero, but there is no evidence for a sequence tending to zero. 
This means that there will be some Dirichlet kernels that seem very concentrated 
near the diagonal (those with tat+i close to zero), but it is unlikely that we can 
improve this behavior indefinitely. 

and [22| illustrate Dirichlet kernels DN{-^y) for fixed y that are not 



Figs. 21 



concentrated near x = y^ while Figs. [23| and [24| illustrate Dirichlet kernels that are 
moderately well concentrated. 
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DKgQ(x,(R-1)/R) 

80 I 1 ! f [ f f 1 f 1 1 

?0 - 

60 - ' 

m - 

40 - 

30 - 

m - 

10 - 
-K - 

-20 1 ^ u \ \ \ \ \ 1 1 

0: &:2 0.3 0.4 0.5 0.6 0.7 0-8 tJ.9 t 

X, Distorted Cantor set 

Figure 21. Dirichlet kernel for CLP with n = 60 and R = 8 
centered at the left most point in the right half of the Cantor set. 
We see that the graph is reasonably small away from the center. 

DK3g(K,(R-1)/R) 

140| \ i \ 1 ■! 1 , 1 , 1 

120' 

80 - 
60 - 
40' 
■20 - 

0.= : : _ - ~™ ^ 

-20 - 

:-40.l 1 ^ ' ' ' ' ' j ' : 

QA 0,2 0.3 0.4 0.5 0.6 0.7 0,8, Q,B 1 

X, Distorted Cantor set 

FIGURE 22. Dirichlet kernel for CLP with n = 89 and = 8 
centered at the left most point in the right half of the Cantor set. 
Again, the kernel is small away from its center. 
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DK5q(x,(R-1)/R) 

70 I , , 1 1 1 , , 1 1 1 



X, Distorted Cantor set 



Figure 23. Dirichlet kernel for CLP with n = 50 and R = 8 
centered at the left most point in the right half of the Cantor set. 
This kernel takes much larger values away from the center than 
the n = 60 case in Fig. [2Tj 



DK^,(x,(R-1)/R) 




X, Distorted Cantor set 



Figure 24. Dirichlet kernel for CLP with n = 99 and R = 8 
centered at the left most point in the right half of the Cantor set. 
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5. Approximate Equalities for CLP Polynomials 



In this section we discuss two types of approximate equahties in the CLP case. 
First we note that by symmetry, Pn{x) is even under the reflection x 1 — x when 
n is even and odd when n is odd. Nevertheless, the plots of P2n{x) and P2n+i(^) 
appear very similar (see n = 98, 99 in Fig. 20). This is especially true on the right 
half of the Cantor set. In Fig. 25 we graph P^oix) and P5i{x) for R = 8. The 
qualitative similarity is striking, but it is difficult to quantify. Fig. [26] shows the 
graph of the difference, and Fig. [27| shows the graph of the ratio. We can give a 
rough explanation using Eq. Q, which for even values may be written 

-1/2 



2n+l 



{x) 



r2n 



2n{x) 



i{x). 



(23) 



^2n+l ^2n+l 

For large values of we may have r2n close to zero and r2n+i close to 1/2. So the 
second term on the right side of Eq. ([23| is close to zero. But on the right half 
of the Cantor set, x is close to one, so the coefficient {x — l/2)/r2n+i is close to 
one. Therefore, Eq. (23) says that P2n+i(^) ^ P2n{x) on the right half of Cr. By 
the odd-even behavior, we have P2n+i(^) ^ —P2n{x) on the left half of Cr. Since 
we see a qualitative similarity of the plots on all of Cr^ we must attribute these 
approximate equalities to an approximate reflectional symmetry across the x-axis 
of the graphs of all Pn {x) . This is roughly apparent in Fig. [20j but it does not hold 
up to close inspection. In particular, the graphs of P2n+i(^) and P2n+i(l — x) are 
not that close (of course P2n+i(l — x) = — P2n+i(^) exactly). 



0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

X, Distorted Cantor set 



Figure 25. Pso{x) and P8i(x), CLP R = S 

The second type of approximate equality refines the above idea. If we express 
the points of Cr as infinite binary decimals, then x (1 — x) simply interchanges 
all digits in the binary expansion. Let T^(x) denote the map Cr Cr that 
interchanges the first m binary digits, leaving all other digits unchanged. In other 
words, Tm permutes the 2^ Cantor subsets of level m by reversing the order of the 



subsets. In Fig. 28 we show the plots of Pi7{x) and Pi7{Ts{x)) for R = 8. There is 
clearly a strong qualitative fit, but the agreement of the two functions is not very 
close numerically. The same pattern persists for P2m_^i{x) and P2m-^i{Tjn{x)) for 
all values of R. At present we have no explanation for this phenomenon. 
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X, Distorted Cantor set 



Figure 26. Pso{x) - P^i{x), CLP = 8 



0.75 0.8 0.85 

X, Distorted Cantor set 



Figure 27. P^q{x) / P^i{x) , CLP R 

Self Similarity: 0P17 and OP17(T4(x)) (R=8) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Distorted Cantor (R=8) 



Figure 28. Pi7{x) and Pi7{Ts{x)), R = S 
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We saw in Figs. 15 



6. CLP Polynomials on Gaps 
and 



17 



that the graph of P2n(^) on the gaps in Cr is 
approximately Gaussian, for n large. In reality, we find that P2n\[i/RXR-i)/R]{x) = 
C2n exp (— (i2n^^^^^) where a{x) = a2n{x) = 2 at x = 1/2, and a{x) ^ 2 for x 
around 1/2. Fig. [29] shows a{x) = a2n{x) for n = 50, CLP R = 4. Due to 
symmetry, we only view a on the interval [1/i?, 1/2]. From the data, we conjecture 
that a2n converges as n ^ oo on any closed interval of the form [1/i? + e, 1/2] to a 



function resembling that in Fig. 29 




Figure 29. a{x) for CLP R = 4, 



50 



The data also indicate exponential growth (in n) of P2n{x) on the central gap 
[1/i?, {R — 1)/^]- This exponential growth results from the behavior of Eq. Q at 
X = 1/2 (with An = 1/2 for all n) as follows: 



P2„(l/2) = (-irn(^2i-i/r2i). 



(24) 



Numerically, we find that {r2j-i/r2j) is bounded above 1 with few exceptions. 
In the CLP R = 2.5 case, there are 134 instances that (r2j-i/r2j) < 1 for j = 
1, . . . , 5000, in the R = 3 case there is one, and for the R = A case there are none. 
Thus, as R increases, {r2j-i/r2j) is bounded above 1 more consistently, as we have 
seen in the R = 8 case in Figs, [g] and [t] Eq. ( [24| ) therefore (generally) gives 
exponential growth in n for |P2n(l/2)|. 

We then summarize the behavior of P2n{x) on the gap [1/R^{R — l)/R] as 
follows: C2n = ^2n(l/2) grows exponentially (according to Eq. |24|, d2n grows 
linearly, and Q^2n(^) appears to converge (to something resembling Fig. [29| ). In 
fact, we can say more: for k G Z+, write 



J 

E 

fc=0 



aj^k{x-l/2Y 



(25) 



Eq. 
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gives all aj^o (observe that a2n,o = C2n and a2n+i,o = for all n), and then 
we can again use Eq. H to compute all of the aj^k- (Here we read Eq. [4] as an 
equation of polynomials?) If we view the coefficients {aj^k} as a (lower triangular) 
matrix, then we see that the left column {o^j^o} and the diagonal {oijj} determine 
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all of the other a^^/c, via Eq. [4] Thus, the behavior of the left column and diagonal 
(and the rj coefficients) dictates the behavior of the other aj^k- 



The left column and diagonal grow log-linearly, as we see in Figs. [3Q|and[3T]and 
Table [l] Note the similarity in the three error plots of Figs. 30 and [31] If we label 
the (top) even indexed error ^evem the (bottom) odd error ^oddi and the (even) 
error function from Fig. [30]^even, we have the following approximate equalities for 
= 0, ...,50: <l>e^en(2fe) ~ <l>odd (2fc + 1) + 1 .36 ^ ^e^en(2/c) + .67. Also, from Table 
[l]we see that sign(aj^/c) = — sign(aj^/c+2), for appropriate j, /c. For the column and 
the diagonal of the matrix of o^j^ki the slope of each best fit line varies with logi?, 
as we see in Fig. [32] 



index] (even) 



Data minus linear fit 



20 



index] (even) 



Figure 30. Logarithmically scaled left column {Q^j,o} of polyno- 
mial coefficients, with error, CLP R = 8. Here we only plot j even 
since for odd j, aj^o = 0. 



a.., CLPR8, slope=2.1863 




Data minus linear fit, even] 



20 40 eo 

index] 







20 40 60 30 100 

index] (even) 
Data minus linear fit, odd] 


V, 





80 100 



index] (odd) 



Figure 31. Logarithmically scaled diagonal {<^j,j} of polynomial 
coefficients, with error, CLP R = 8. 



On the other gaps (i.e. the images of [l/R^ {R — 1)/^] under compositions of 
the maps Fi and F2), we observe similar qualitative behavior. That is, on these 
smaller gaps, Pn{x) resembles a Gaussian for all large and even n. We have seen 



this self-similar behavior already in Figs. 15 and 16 However, the local extremum 
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Slope of best fit line of (log(a. .) vs. j) , and for (log(a. ^) vs. j (even)), for CLPR 



log(R) 



log(R) 



Figure 32. Slopes of best fit lines for ajj and aj^o (as in Figs 



30 and 31 ), for G {2.5, 3, 4, 8, 16, 32}. 
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-26134 





26237 




Table 1. 
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= 8, Matrix of (rounded) aj^k 


coefficients, di 









vided by 10^ 



Local Extremum of P (x) on [1/R^(R-1)/R^], CLP R=4 



Figure 33. x-value of local extremum of Pn{x) on [1/R'^^{R 
l)/i?2] = [1/16,3/16] as a function of n, CLP R = 4. 



of P2n{x) on [1/i?^, (i? — I) / R'^] is uot in general the center of the given gap. We 
glimpse this phenomenon in Fig. 



17 In Fig. 



• 2n|[l/16,3/16] 



in the CLP R = 4 case. The center of this gap occurs at x 



33| we plot the local extremum of 

.125, 



but the local extremum fluctuates with average around .122. 
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7. WLP AS A FUNCTION OF n FOR FIXED X 

For the classical Legendre polynomials on [0,1], we find structure in the val- 
ues of Pn{x) as functions of n for fixed x. As an example, for x = 0, Pn{x) = 
(-l)^V2n + 1, and for x = 1, Pn(l) = V2n + 1. The computations of Pn(0) and 



of Pn(l) are shown at the top of Fig. 34 These algebraic relations are a renormal- 



ization of the behavior of the classical Legendre polynomials on [—1, 1], where 
werecahthat Pn(l) = 1, Pn(-l) = (-1)'' and | |Pn| Ii2[_i^i] = l/(n + l/2) forn > 1 
IGautschi 20041 . We now increase the WLP weight pi at the points x = 0, 1. As a 
result, we see a perturbation of the behavior of the classical Legendre polynomials 



(on [0, 1]) in Fig. 34 For pi .7, a transition seems to occur in the behavior of 
{Pn(0)}^x. At this transition, we appear to have a multiplicative periodic function 
Pn(0), i.e. a function / : Z+ R where for some a > we have f{ax) = f{x) for all 
X. For pi = .7, we can logarithmically scale the x-axis as in Fig. [35] to see a nearly 
multiplicative periodic function. For other pi, we evidently have a multiplicative 
periodic function P^(0)/n^ for an appropriate choice of P = P{pi). 

For generi([^ X G [0,1], if we plot the vectors (Pn(x), Pn+i(x)) in the WLP 
Pi = .5 (i.e. classical) case, then these vectors are attracted to an ellipse as n ^ 
oo. This ellipse is centered at the origin, and its axes and orientation vary with 
X. Also, the vectors (Pn(x), Pn+i(x)) rotate around the ellipse in a periodic way. 
That is, for a certain k > 1 which depends on x (where /c G Z is the period), 
(i((Pn(x), Pn+i(x)), (Pn+/c(x), Pn+fe+i(^))) < e for some small e > 0. Geometrically 
speaking, (Pn(x), Pn+i(x), n) G travels along a helix near the surface of an 
ellipsoidal cylinder. Thus, plotting (P^(x), Pn+i(x)) projects this helix onto the 
plane. In Section [sj we focus on the behavior of {Pn{x)^ Pn-\-i{x)) for generic x G 
[0,1]. For now, we merely note the contrast between the images of generic and 
no n- generic points. 

As observed for x G {{0}, {1}} (a set of two non-generic points), increasing pi 
perturbs the dynamics of the vectors (Pn(x), Pn+i(x)). Therefore, increasing pi 
should result in perturbed dynamics for other non-generic x. To this end, consider 
X G {{.25}, {.5}, {.75}}. For the classical Legendre polynomials on [0,1], these 
X values yield finite attractors for {(P^(x), Pn+i(x))}^;L- These attractors have 
three, four and six points, respectively. We call the corresponding number of points 
kx- The vectors (Pn(x), Pn+i(x)) travel in a clockwise fashion about the kx attract- 
ing points, as n increases. When we increase pi, the dynamics of {Pn{x) ^ Pn-\-i{x)) 



for X G {{.25}, {.5}, {.75}} change dramatically, as we see in Figs. 36 and 37 In- 



stead of cycling ever closer to kx attracting points, the vectors (Pn(x), Pn+i(x)) 



cycle through kx fractal spiral arms, as we see in Fig. [37| However, in the case that 
pi > .5, the attractor of these spiral arms is difficult to determine. By measuring 
the distance of {Pn{x)^ Pn-\-i{x)) from (0,0), it seems that the spirals are bounded 
away from the origin as n ^ oo. Therefore, by the symmetry apparent from Figs. 
[36|and[37l (Pn(x), Pn+i(x)) should be attracted to a (fractal) ellipse, or a finite set 



of kx points. As a further contrast to the pi = .5 behavior, we see from Fig. [37 
that (Pn(x), Pn+i(x), n) G travels along a "fractal helix" with radius decreasing 
in n. 



By a generic point we mean a point chosen, with respect to a uniform distribution, among 
a suitable set of (rational) values in floating point arithmetic. In all cases, we only treat x 
representable in double precision floating point arithmetic. 
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Figure 34. Pn(x) values vs. n for x = 0, 1 plotted for the WLP 
Pi = .5, .6, .7, .8, .9 families 
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p (X) values for x=0 



log( Iteration number n) 



Figure 35. Pn(0) values vs. n logarithmically scaled in the x-axis 
for WLP pi = .7 
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Figure 36. Pn{x) values vs. n for x = .75, .5, .25 plotted for the 
WLP pi = .6, .7, .8 families. Note that for pi = .5, we would get 
solid horizontal lines. 
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Figure 37. Vectors {Pn{x), Pn^i{x)) for x = .75, .5, .25, plotted 
for the WLP pi = .6, .7, .8 families 
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8. WLP AND CLP Dynamics 
As promised, we now examine the dynamics of Pn{x) for generic x. Here we 



examine both WLP and CLP, and we summarize the results in Figure 38 For 
WLP, we recah for pi = .5 and fixed, generic x G [0,1], there exists an integer 
kx > I so that 

d{{Pn{x),Pn+l{x)),{Pn+kA^),Pn+k^ + lix))) < 6 (26) 

for some smah e > 0. Here we mean that the sequence of points 
{{Pn-^jk^{x)^ Pn-\-jkoo-\-i{^))}j^o is periodic of period up to a smah error of e 
at each step. We therefore make kx the minimal positive integer satisfying our 
condition [26| If we plot Pn{x) vs. n, we find that Pn{x) is a superposition of 
kx sinusoidal functions (see entry (1,1) in Fig. [38|. These kx (approximately) 
periodic functions are given by P(/c^n+j)(^) foi" j ^ {1, 2, . . . , /c^;}. For example, 
we can see that /c.95 = 14 by counting the number of distinct periodic functions 
in the plot of Pn{x) versus n. Since we essentially have a superposition of phase 
shifted cosines, it follows that the distribution of values of Pn{x) is the function 



1/y^1 — suitably rescaled (see entry (2,1) in Fig. 38). Finally, if we color the 
iterates (Pn(x), Pn+i(x)) in a 14-periodic manner, we get entry (3,1) in Fig. 38 
In this plot, the point (Pn+i(^), ^n+2(^)) overlaps all points of lower index. 

So, what happens when we increase pi? As in Section [7| we observe a per- 
turbation of the pi = .5 behavior. This perturbation becomes more exaggerated 
the larger pi becomes. In the middle column of Fig. |38j we show the same plots 
described in the previous paragraph, for = .6. We again choose x = .95 so 
/c.95 = 14, and we observe fourteen vaguely periodic functions. In other words, for 
jf G {1, 2, . . . , kx}^ we see that P(/g^^+j) (.95) is a perturbed sinusoidal function of n. 
We can see these perturbations via the distribution function (entry (2,2)) and the 



plot of (Pn(x), Pn+i(x)) (entry (3,2)) in Fig. 38 That is, the distribution of values 
of Pn{x) is a perturbed version of the function 1/ ^^1 — ?/^, suitably rescaled. Also, 
the attractor of the vectors (Pn(x), Pn+i(x)) appears to be a thickened (fractal) 
ellipse. 

The behavior for generic x in the CLP case is displayed in Fig. [38] in the right- 
most column. As an example, we view the iterates of x = 1/5 = l/{R-\- 1) G Cr 
for the case R = A. (Recall that Pn{x) grows exponentially for x G [0,1] \ Cr^ 
so we expect more information by viewing x G Cr.) For generic x, the distri- 
bution of Pn{x) values is unclear (Fig. |38j entry (2,3)). Nevertheless, the Pn{x) 
values do seem concentrated around the origin. Also, the vectors (Pn(x), Pn+i(x)) 
obey a 2-periodic behavior. Specifically, the points {(P2n(^)7 ^2n+i(^))}^i clus- 
ter along the line {{x^y) G M? : y = —x} (to abuse notation). Also, the points 
{(P2n+i(^), P2n+2(^))}^i clustcr looscly arouud the origin. To observe this be- 



havior, note entry (3,3) in Fig. 38 



In our investigation, we have found only one non-generic point x G Cr H [0, 1/2] 
for the CLP case (other than the endpoint x = 0). (Note that the symmetry of 
these Pn{x) across x — 1/2 reduces our investigation to Cr H [0,1/2].) Unlike 
the generic iterates which demonstrate only 2-periodic behavior, the non-generic 
point X = 1/R exhibits 4-periodic behavior. The columns of Fig. [39| correspond 
to different pi values and the iterates oi x = 1/R. For x = l/R, the vectors 
{P4n-\-j{x)^ P4n-\-j-\-i{x)) form four disjoint attractors for j = {1,2,3,4}. The set 
{{P4n-\-j{x), P4n-\-j-\-i{x))}^=o foi" J = 1 rcsidcs in the lower left quadrant, for j = 2 
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Figure 38. Summary of Generic point values: Pn values, their 
histograms, and the vectors {pn{x)^Pn-\-i{x)) plotted for the Le- 
gendre, WLP pi = .6 and CLP, R = 4 polynomial families 

this set lives in the upper left quadrant, for j = 3 the upper right, and for j = 4: 
the lower right. 
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Figure 39. Pn values, their histograms, and the vectors 
(Pn(x),Pn+i(x)) plotted for x = 1/R, CLP R = 2.5,3,4 and 8 
polynomial families. 
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9. Concluding Discussion 

Experimental mathematics allows us to see patterns that otherwise would never 
be observed. At times these patterns allow us to formula explicit conjectures, 
and eventually proofs may be found to turn conjectures into theorems. But often 
the story is more complicated. The observed patterns may turn out, on closer 
inspection, to be only approximately true. No clear cut conjectures emerge that 
are truly supported by the experimental evidence. In such cases, what should we 
do? 

We could simply discard the experimental results and move on. But the approx- 
imate patterns that we observe may be of great interest. The experiment might be 
trying to tell us something that we are not quite able to capture in the conventional 
format of mathematical statements. Experimental science is often messy in exactly 
this way, and we think it would be a shame to limit experimental mathematics to 
only the clean paradigm. 

We see the results reported in this paper in exactly this light. We do have clean 
cut conjectures on the values in Section |2] and boundedness for WLP in Section 
[3j but beyond that we only have messy evidence. We see some CLP Dirichlet ker- 
nels that look more like approximate identities than most Dirichlet kernels. We are 
able to relate this to the small size of the coefficients. The experimental evidence 
suggests that we will not find a sequence with converging to zero. So we cannot 
offer a conjecture analogous to the results of [Strichartz 2006] about uniform con- 
vergence of certain special partial sums of the orthogonal polynomial expansion of 
an arbitrary continuous function. However, the results of [Strichartz 2006 suggest 
that perhaps such a statement might hold only for special values of and we have 
only examined a few choices of R. 

The approximate equalities for CLP discussed in Section |5] appear very striking 
when one looks at the graphs of the polynomials. It is on closer inspection that 
one observes the deviation from equality. A clean conjecture might say that the 
deviation goes to zero in some limit, but the evidence does not really support 
such a conjecture. Nevertheless, we find the approximate equalities intrinsically 
interesting. 

Similarly, a glance at the graphs of CLP on the gaps suggests they are close 
to being Gaussian. Of course one comes to expect Gaussian limits in many math- 
ematical situations. But here in Section |6| on closer inspection, we see a decided 
deviation from the Gaussian model. In this case we believe there is a true limit 
law; we just have not been able to find it. 

Perhaps the most interesting contribution of this paper is the dynamical per- 
spective discussed in Sections [7| and [S] Here we may be excused from offering 
explicit conjectures on the grounds that the pictures offer a view of significant 
structures that were previously unrecognized. We hope that others will be inspired 
to investigate the dynamics perspective in a systematic way. 
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